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Abstract. Isospin-breaking effects in nuclear matter are studied in the framework of a medium-modified 
Skyrme model. The proposed effective Lagrangian incorporates both the medium influence of the sur- 
rounding nuclear environment on the single nucleon properties and an explicit isospin-breaking efi'ect in 
the mesonic sector. The approach predicts that the neutron-proton mass difference decreases in isospin- 
symmetric nuclear matter but by a very small amount only. 
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. 1 Introduction 

■ The evaluation of isospin-breaking effects in a dense nu- 
' clear medium is an interesting problem in nuclear physics. 
. This is particularly the case for the Nolen-Schiffer anoma- 
' ly [I] in mirror nuclei which may be explained by a change 
, of the neutron-proton mass difference in the nuclear en- 

■ vironment. Henley and Krein [2] argued that the Nolen- 
\ Schiffer anomaly can be resolved if the neutron-proton 
r mass difference decreases rapidly in the medium. For this 

reason various approaches for calculating the neutron-pro- 
ton mass difference in a dense nuclear environment have 
' been proposed in the hterature [3l|4l[5ll6l[7ll8l[9] . Some of 
these approaches indeed confirmed a decrease [3l|4l|5l[6l 
[7] of the neutron-proton mass difference in the nuclear 
medium, while others gave opposite results [8l[9]. Since 
the Skyrme model [lOlfTT] treats the properties p ^ll3U14j 
and the interactions [T51ll6l[T71[T5] of the nucleons on an 
equal footing, it is interesting to study what this model 
has to say for possible medium-modifications of this quan- 
tity. Note also that the neutron-proton mass difference is 
composed of a strong and an electromagnetic contribution 
of comparable size, which might be affected differently in 
a dense medium. 

In the present work we evaluate isospin-breaking ef- 
fects in the baryonic sector of the Skyrme model by incor- 
porating the influence of the surrounding nuclear environ- 
ment on the nucleons as e.g. in Refs. p!9ll20j . In order to 
perform our studies we modify the nonlinear cr-model La- 
grangian, which is generalized by explicit isospin-breaking 
effects in the mesonic sector [21] , based on the well-known 
pionic field equations in the nuclear medium [22j . Further- 
more, the usual Skyrme term is added in order to stabilize 
the solitons of the model. 



The paper is organized as follows. In Sec. [2] we discuss 
our model Lagrangian. Section [3] is devoted to the study 
of the field equations to be solved. After an explanation 
of the quantization procedure and the derivation of the 
strong part of the density-dependent neutron-proton mass 
difference in Sec. HI we will present the electromagnetic 
(EM) part of that mass difference in terms of electromag- 
netic form factors in Sec. [5] The input values of the model 
parameters used in this work are discussed in Sec. [H Sec- 
tion [7| is devoted to the presentation and discussion of the 
results. Finally, in Sec. [5] our conclusions arc summarized 
and an outlook for further studies is given. 



2 Medium-modified Lagrangian 

We start from a generalized Lagrangian which incorpo- 
rates an explicit isospin breaking term in the mesonic sec- 
tor 

C = C2 + Ci^ + /^gxSB , (1) 

/:2 --^Tri,,L^ (2) 
A = ^Tr[i^,L,]2, (3) 

/:gxSB = {Tr {U - l)M+{U^ - l)M + 

-Tr [U - l)M-T^{U^ - l)M-T^] , (4) 

where = U'^d^U . The Einstein summation conven- 
tion for repeated Greek or Latin indices is adopted, unless 
stated otherwise. £2 is the usual Lagrangian of the nonlin- 
ear sigma model, £4 is Skyrme's stabilizing fourth order 
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term. Cg^SB is a generalized pion-mass term proposed by 
Rathske 21 which includes explicitly isospin-breaking ef- 
fects in the second term. Note that the latter has the same 
structure as the strong isospin-breaking fourth-order La- 
grangian of chiral perturbation theory, see the term pro- 
portional to the low-energy constant L7, Eq. (6.16) of 
Ref. 23J. El Skyrme models with isospin-breaking terms 
proportional to the p-to meson-mixing, as e.g. discussed in 
Refs. [24ll25j . are actually of the same order in the chiral 
expansion, but more involved because of the additional 
meson degrees of freedom. The parameter F^^ is the pion 
decay constant, while e is the dimensionless Skyrme con- 
stant. The SU(2) matrix U{t,r) = exp{2ir • ^{t,r)/F^} 
represents the isotriplet of pion fields Their masses, 
7Tijr± and m^o, enter the Lagrangian in the following com- 
binations 

M± = ^iml^±ml„)/2. (5) 

The medium-modified Lagrangian, which takes into 
account the influence of the surrounding nuclear environ- 
ment on the pion fields, can be set up by following the 
steps presented in Ref. [19]. It has the forn:|3 

£* = £2 + A + ^2,gxSB , (6) 

^2,gxSB = -§{ - (X^ + Xl) Tr VU ■ VC/t 
+ (Xp -Xp) TrVf/ra- VC/trg 
+ {2Ml + xt + Xs) Tr(C/ - 1){U^ - 1) 
- {2Ml^ + xt - Xs) Tr(t/ - l)r3(C/t - l)r,}. 

Here, xtp Xsp functionals of the low-energy S- 
and P-wave pion-nucleon scattering lengths and volumes, 
60,1 and co,i, respectively [22]. They express the influence 
of the medium on the charged and neutral pion fields, 
TT^ and tt". Furthermore, they depend on the densities of 
the neutron and proton distributions, and pp, of the 
environment. One can check that the Lagrangian ([S]) - in 
the linear approximation - generates the well-known field 
equations (9^9'" -|- m^^ „ + 77±'°)7r=^'° = [22], where the 
polarization operators IJ^'^ have the schematic form 

iI^'° = Xs^^° + Vx^^°-V. (7) 

The medium-modified Lagrangian ([6]) can be rewritten 
in the following way 

= £4 + /:*SB + > (8) 
= ^{TraoC/aoC/t-(l-x°)TrV(7- VC/t}, (9) 

^ This is manifest when the isospin-breaking term is rewrit- 
ten as in Eq. Because of UU^ = 1, Rathske's first term, 
the pion-mass term, is compatible with the usual choice of ex- 
plicit chiral symmetry breaking proportional to Tr {U + U^ —2) 

m- 

^ From now on an asterix indicates a density-dependent 
quantity. 



'CxSB = -^^(l + ;^jTr(^-l)(t/^-l), (10) 

AH* ^^^[{X^-X\:) [Tr(r.VC/).Tr(r.V[/t) 
-Tr(r3VC/)-Tr(T3Vf/t)] 
- (2M2_ + ^± _ [Tr (t,C/) Tr(T,C/t) 

-Tr(r3[/)Tr(r3[/t)] }, (11) 

where £^gB is the chiral-symmetry-breaking term and AL* 
is the isospin-breaking term, separately. Note that in the 
absence of AL* the Lagrangian ([H]) reduces to the medium- 
modified Lagrangian presented in Ref. [T^] • 

In general, the strong isospin breaking in the bary- 
onic sector is generated by two effects: an explicit isospin- 
breaking in the mesonic sector (in the Lagrangian this 
term is proportional to A^?_) and the isospin-dependent 
influence of the medium on each of the pion fields, Xp ~ 
Xp and x^ ~ Xs - For simplicity and in order to concen- 
trate on the effects due to explicit isospin breaking in the 
mesonic sector, we will consider from now on only isospin- 
symmetric nuclear matter: 

±^ 0= 47rcop/7? 

Xp - Xp ^ ^ 4g'7rcop/r; ' 

Xs"^ = Xs = Xs = -47r?7(6o)offP • (12) 

Here ry = 1 + m^a /m^ is a kinematical factor, = 
938 MeV is the nucleon mass, g' is a correlation parame- 
ter, Co is the P-wave isoscalar pion-nucleon scattering vol- 
ume, (&o)cff is the corresponding effective S-wave isoscalar 
scattering length in nuclear matter, and p is the density of 
the surrounding medium [22]. Moreover, we will limit the 
investigation to homogenous nuclear matter, Pp = Pn = 
p/2 — const, and to the evaluation of the neutron-proton 
mass difference in nuclear matter, Am^^ = rn* — m*. 

To finish this section, we remark that the approach 
used here should be distinguished from the method used 
in Refs. [2511271125] . In these papers, all meson properties 
like the masses, decay constants, etc. were assumed to be 
density-dependent and taken from a microscopic model 
of the meson dynamics. Here, the density-dependence is 
parameterized in terms of phenomenological input, more 
precisely the effective scattering length (&o)eff and volume 
Co- 

The Lagrangian ([6]) is of course only one representative 
of a larger class of Skyrme- type models with e.g. higher- 
order terms or vector-meson degrees of freedom which, 
rooted in a hadronic language, do allow for a simulta- 
neous description of nucleon properties and interactions. 
However, the present choice with its scale-independent 
dimensionless Skyrme parameter e can be considered as 
the simplest, but also as a generic one, with a minimal 
number of derivatives and model-parameters, that incor- 
porates explicitly strong isospin breaking, stable solitons, 
and medium modifications. 
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3 Classical solitonic solutions 

Due to the isospin-breaking term of the Lagrangian 
the third component of the isovector field is singled out in 
the field equations^: 



dC* 



dC* 



0. 



(13) 



Treating the (strong) isospin-breaking as a small perturba- 
tion, we can distinguish between a non-perturbed baryonic 



background with 



and perturbative contribu- 



tions proportional to the very tiny strong-interaction in- 
duced part of m^± — m^o ^ 0. The non-perturbed system 
has the usual time-independent classical (baryon number 
-8 = 1) solutions #(t — 0, r) = (t^if) of the isospin- 
symmetric Lagrangian, whereas a perturbation around the 
classical background <p{r) can be absorbed e.g. by a time- 
dependent modification. 

We follow Rathske [5T] by using for this purpose the 
separation ansatz 



#(t,r) =r(t)0(r) 



(14) 



where a manifest isospin-breaking (with respect to the 3rd 
isospin-axis) is encoded in the time-dependent iso-rotation 
matrix T: 



cosa*t sma*t 
T{t) = I -sina*i cosa*t 
1 



(15) 



Here a* is a density-dependent global constaniQ which 
would vanish in the absence of the small perturbation 
ZirriTr ~ m^± — m^o^ since then the system would reduce 
to the non-perturbed configuration lim C* = C^p- 

Anij^—^O 

The ansatz (|14p together with (|15p . when inserted into 
the Lagrangian ([8|) and after spatial integration, allows to 
cancel the non-perturbative shift (IB-NP) proportional to 
TO^± — m^o that results from inserting the non-perturbative 
(NP) background (p{r) of the isospin-symmetric case into 
the isospin-breaking part (IB) of the Lagrange function 
against the isospin-breaking term (IB-T) that results from 
the time-dependent rotation by the T matrix: 



L* = J C* d^r = -A/*p - MlA*_ + —A* 

= I C*^p dV + I A*B-NP + I A*B-T . (16) 

Here M^p[0] is the isospin-symmetric mass of the non- 
perturbed configuration. The moment-of-inertia type quan- 
tities A* = A*[4>] and A*_ = A*_[(p], belonging to the T- 
rotation and the non-perturbative shift, respectively, are 
functionals of the pion fields. The global parameter a* 



^ The explicit form of Eq. (|13p can be found in Ref. [21) for 
the free space case. Note that the Skyrme term is omitted there. 

* Note that a sign change of t can be compensated by a sign 
change of the angular velocity a* . 



serves here as a constraint parameter, whereas in Ref. |21j 

- under the neglect of stabilizing higher-order terms - the 
ansatz (|14p also extremized the remaining action. Thus 
the value of a*^ is fixed here by the following condition 
which implies the cancellation between the IB-NP and IB- 
T parts in the Lagrange function: 

A* 

a*^ = 2M^—. (17) 

^ A* ^ ^ 

As the iso-rotation matrix T is chosen in such a way that 

- at the classical level and after spatial integration - the 
isospin-breaking is rotated out of the system, the clas- 
sical field equations (fT3|l effectively reduce to the non- 
perturbed isospin-symmetric case: 



dC* 



dC* 
'd¥ 







5M*p[cl)] 



0. (18) 



In general, one should consider non-spherical time- 
independent field configurations (p{r) if the density gra- 
dients of the surrounding nuclear environment are large. 
In this case, the static field configuration acquires a 6- 
dependence [SU], i.e. (p — 4>{\r\,9), even if isospin-breaking 
terms are absent. The above presented "rotation" pro- 
cedure has the advantage that it remains applicable for 
very non-uniform density-profile-dependent cases, in other 
words, for finite nuclei. But as stated in Sec. [21 only the 
case of an isospin-symmetric homogeneous nuclear envi- 
ronment with a constant density is considered here. Thus 
the use of a spherically symmetric static configuration 
4>{r) — {r/r)FT^F{r)/2 is still appropriate. Under this 
hedgehog ansatz the terms defined in (I16|) are given by 



M^p = TT 



Fl 



f(i-Xp)M^, 



2S^ 



+F^ (mlo + Xs) (1 - cos^^) dr , (19) 

OQ 

^* = |/k + lfe + 5)U^r^d.,(20) 



(21) 



where F^ = dF(r)/dr. Also the abbreviation S = sinF(7') 
is used here. 

The corresponding field equations (fTS]) do not explic- 
itly depend on the isospin index and therefore reduce to 
the simple radial form that is used in this work to deter- 
mine the profile function F{r): 

Fl , f 2 sin2F 



2 (2w^F 



sin2F / , sin^i^ 
Frr + ^ Fl ^ 



F, 



I (m^o + X.) sinF = . 



(22) 
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Here Frr stands for the second derivative of the profile 
function F(r) with respect to the radial coordinate r. 

If the Skyrme term were omitted, the condition (jl7p 
would reduce to a*^ = w? ± — mio and the results of the 
work [21j would be reproduced in a natural way. However 
in this case, due to the then profile-independent value of 
the parameter a*, the density dependence of the isospin- 
breaking effects would be losll£|. From this point of view, 
the presence of a higher-order term in the Lagrangian 
is crucial, not only for its stabilization role, but also for 
the description of the density- dependence of the isospin- 
breaking in the nuclear environment. Of course, there is 
a model-dependence due to the nature of the stabiliz- 
ing term. We reiterate that the Skyrme term is the most 
simplest choice, with the least number of derivatives and 
model-constants, that guarantees the existence of stable 
baryons. In fact, the chiral order of the Skyrme term is 
the same as the one of the strong isospin-breaking term of 
chiral perturbation theory. 

To finish this section, let us state the expectation that 
the above mentioned shift induced by the isospin-breaking 
terms and the compensating isospin-breaking rotation will 
be manifest in the quantized theory, although they are 
constrained to cancel each other at the classical level. 



4 Quantization and the strong part of Am 



np 



exp i— T • ^{t, r) 



The standard quantization procedure of the Skyrme mo- 
del [12] requires time-dependent rotations in space, TZ{t), 
and isospin-space, T{t): 



U{t,r) 



^U'{t,r')^expi^i—T-I{t)^{t,n-\t)r)^ ,(23) 

where r) is defined in Eq. ^1)1 . The rotations TZ{t) 
and iso-rotations X{t) correspond to standard collective 
zero-energy modes of the classical soliton, whereas the ma- 
trix T of Eq. describes a constrained rotational mode 
with respect to the classical soliton. The SO (3) matrices, 
I{t) and Tl{t), satisfy the conditions 

(t) = ^^3l^l , T^7k (^)^fcj (*) = -^vi ' (24) 

where cui and f2i are the angular velocities of the isospin 
and spatial rotations, respectively, and the 'dot' symbol- 
izes a time derivative. Under these rotations and the con- 
straint (fT7|) the corresponding Lagrange function L* = 
J C* d^r takes the form 



-Af, 



NP ■ 



A* 

T 



{uj - fiy + 2a* (w3 - J73) . (25) 



Note that the third components of the angular velocities 
couple to the constrained angular velocity a* . Finally, the 



^ The isospin-breaking effects resulting from the pion-mass 
variation in the nuclear medium are beyond the scope of the 
present model. 



usual Legendre transformation utilizing the canonical-con- 
jugated operators Ti = duj^L*, Ji — dn-L* of LOi and Qi, 
respectively, leads to the quantum Hamiltonian: 



H* 



M, 



NP 



NP 



2/1* 



MlA* 



T2_ 

2A* 



T 

2A* 



Ts - A*a' 
2A* 

- a*T3 . 



(26) 



There does not exist a canonical-conjugated momentum of 
the angular velocity a* since this is a constrained quantity. 
The mass-shift term due to the explicit isospin-breaking, 
the second term of Eq. (|26p . reappears - compared with 
Eq. (fTO|) or Eq. (|25p - after applying the constraint pT|) 
to the induced term {A*/2)a*'^. Moreover, note that the 
isospin breaking in the mesonic sector is manifest at the 
quantum level by the coefficient of the third component of 
the isospin operator, see the fourth term of Eq. (|26p. Con- 
sequently, by evaluating the quantum Hamiltonian (|26p 
between appropriate nucleon states one can isolate the 
strong part of the neutron-proton mass difference in the 
nuclear medium as|f| 



*;(strc 



^*(strong) ^ /ij„*(stro„g) ^ ^* 



(27) 



This term is enhanced in comparison to the tiny mass shift 
M^A'^, since the latter is quadratic and not linear in the 
small parameter a*. 

Finally, let us discuss the scaling under the expansion 
in the number of colors Nc of QCD, the underlying the- 
ory of strong interactions. From the viewpoint of the non- 
relativistic quark model, the mass difference of a large 
Nc neutron (with (Nc + l)/2 d-quarks and {N^ - l)/2 
w-quarks) and a large Nc proton (with {Nc + l)/2 u- 
quarks and {Nc — l)/2 d-quarks) is expected to scale as 
TTT-d — iTiu ~ Nc - Let us compare this with our Eq. (|27|) : 
As usual, we start with the assumption that Ft^ ~ 1/e ~ 
\//Vc, such that the soliton mass Mj^p, the moment-of- 
inertia A* and also A*_ scale as 0{N^). Assuming that 
Ai- scales as A^+, namely as 0{Nc), Q the constrained 
constant a* (sec Eq. ([T7|) and the strong neutron-proton 
mass difference ^7} is of the same order 0{N^). There- 
fore, the Nc scaling of Eq. (pS)) is as it should be: the 
non-perturbative isospin-symmetric mass term M^p and 
the non-perturbative (via isospin-breaking induced) mass 
shift A4^A*^ scale both as 0{N^); the quantization term 

T /(2yl*) scales as 0{Nc^), as expected for the quanti- 
zation of zero-modes of a soliton [12]; finally the isospin- 



^ Note that the sign of a* is now fixed by the known sign of 
the strong nucleon-proton mass difference in free space. 

^ Note that the Nc scaling of the SU(3) low-energy constant 
L7 ~ (l/48)F^/m^, ~ 0{N^) [23] would actually imply that 
should scale as 0{Nc ) (see ((4|), since the rj' would become 
a Goldstone boson in the limit Nc ^ 00 {i.e., the rj' mass 



scales as m„ 



For Nc = 3, however, the prefactor 



is very tiny, such that the above assumed 0{Nc) behavior of 
A^i = {m'^± — m^o)/2 is the more natural choice for the strong 
pion mass difference. 



Ulf-G. Meifiner et al.: Neutron-proton mass difference in nuclear matter 



5 



breaking term a*T3 resulting from the constrained iso- 
rotation (with the enhancement by A* ^ Nc relative to 
the quantized isospin, see the first hne of Eq. ([26])) scales 
as 0{N^). 

5 Electromagnetic form factors and the 
electromagnetic part of Am*p 

The electric (E) and magnetic (M) form factors of the 
nucleon are defined through the expressions 

Gm(9') = 5 ™n 1 dV e^-^-ir x j{r)] , (28) 

where is the squared momentum transfer. Furthermore, 
and j correspond to the time and space components of 
the properly normalized sum of the baryonic current B* 
and the third component of the isovector current of 
the Skyrme model, i.e. 

S; = ^£M.a/3TrL^i"L^ (29) 

V^''* - ^ Tr jra [-F^ C,L, + 1 [L\ [L„ L.]]) | 
+ {L^R), (30) 

c-={i4:;::;:2.3. m 

Here e^^ap is the totally antisymmetric tensor in four di- 
mensions and the definition i?^ = Ud^U'^ is used. Note 
that the index /i in Eq. (j30p is not summed over. By evalu- 
ating these current operators between appropriate nucleon 
states, one obtains the density-dependent electromagnetic 
form factors of the nucleon p6l[271 28,29, 30 . 

Although the isoscalar (S) electric and magnetic form 
factors of the nucleon 

oo 

Gl* ^-^jPr sm^Fjoiqr)dr, 



oc 

^s* "IN /■ J, . 2 p 2 jijQr) , . 

Gm = / sm 1< r dr (32) 

ttA* J qr ^ ' 



do not explicitly depend on the density, there still exists 
an implicit medium modification caused by the density 
dependence of the profile function F. This holds also for 
the isovector (V) electric form factor 

oo 


'K sm^ F jo{cir) dr . (33) 



The isovector magnetic form factor, however, has both 
explicit and implicit density dependences: 

oo 


xsin^i^r^^iMdr. (34) 
qr 

The medium-dependent form factors of the proton and 

neutron are defined as G^{^ — G*!*,^ ± G^*^j with the 
normalization conditions &^ [0) = 1^ Gg*(0) = 0, GPi*(0) 
= /ip, G^(0) = /i*, where /i* and /i* are the in- medium 
magnetic moments of the proton and neutron, respec- 
tively. Note that the magnetic moments are measured in 
terms of the nuclear Bohr magneton (abbreviated as n.m. 
in Tab.[T]) with the free-space nucleon mass m^. Therefore, 
Ton and not to^ is used in Eqs. ([15]), and ((M)) . 

Finally, for calculating the electromagnetic part of the 
neutron-proton mass difference, one can apply the formula 

oo 


see e.g. Ref. [3T]. Here a = e^/47r « 1/137 is the elec- 
tromagnetic fine structure constant (with e the elemen- 
tary charge). Note also that the right-hand-side of this 
equation, after Eqs. ([5^ to are inserted, does not de- 
pend on the value of the nucleon mass and moreover that 
Am*^^^^ scales as 0{N^,) hke Z\TO*[f'""s). 



6 Input parameters of the model 

As input for the free mass of the neutral pion we take the 
PDG- value [32]: to^o = 134.977 MeV. This choice induces 
the values = 108.11 MeV and e = 4.826, if one insists 
on reproducing the empirical (isospin-averaged) masses 
of the nucleon and delta, ton = 938 MeV and Ma = 
1232 MeV, in free space {p — 0) and without isospin break- 
ing term (A4_ = 0). The small differences to the values 
F„ = 108 MeV and e = 4.84 of Ref. [l3| are caused by the 
value TO^n-o = 138 MeV used there. 

The mass to^± of the charged pions, actually only its 
strong contributiorQ to AttIt^ = my^± — m„o , is extracted 
as a variational parameter from the fit to the empirical 
value Z\TOnp'^'''' = 1.29 MeV in free space. To do this, we 
first calculate the electromagnetic part of the neutron- 
proton mass difference in free space in analogy to the 
formula dSS]), Z\mi™^ = -0.68 MeV. This is in good 
agreement with the result of Ebrahim and Savci [33j ob- 
tained long ago. Model calculations, which include tt, p 

* The electromagnetic contribution to m^± is not considered 
throughout this paper. 
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Table 1. Masses and magnetic moments of the nucleons in 
nuclear matter with a density p = Apo, where po = O.Sm^o is 
the saturation density of ordinary nuclear matter. 



A 


ml [MeV] 


ml [MeV] 


Mp [n-m.] 


Hl [n.m.] 





937.4 


938.7 


1.99 


-1.27 


0.2 


861.9 


863.2 


1.86 


-1.13 


0.4 


796.6 


797.9 


1.75 


-1.02 


0.6 


739.6 


740.9 


1.65 


-0.92 


0.8 


689.5 


690.8 


1.56 


-0.83 


1.0 


645.2 


646.4 


1.48 


-0.75 



and Ld mesons, also lead to similar results [S^. Moreover, 
these results also are consistent with the value zimi™-' = 
(—0.76 ± 0.30) MeV estimated from the Cottingham for- 
mula in Ref. [31j . Finally, applying 



(stron. 



(36) 



we find the following estimate of the neutron-proton mass 



difference due to the strong interactions 

,(Exp) 



(strong) 



np — 

^TOnp"''^ - ^m'np'' = 1.97 McV which is fitted in the 
present model by adjusting (the strong part of) AttIt^, i.e. 
the quantities A4+ and Ai- oi Eq. ([S]). In contrast to the 

sizable value of Z\mnp'^°"^^ the pion-mass difference due to 

strong interactions is very smalQ Z\mt*'*™"®'' ~ 0.04 MeV, 
compatible with the results of Ref. [31] . 

Finally, the parameters of the pion self-energy are taken 
from Ref. [22], Table 6.2: g'^ = 1/3, cq = 0.21m-^ and 
('>o)eff = — 0.024m~^. The latter two are based on the 
tabulated isospin 1/2 and 3/2 S1/2, P1/2 and P3/2 pion- 
nucleon scattering lengths from Ref. [35j which are still 
compatible with modern calculations in the framework of 
chiral perturbation theory, see e.g. Tab. 2 of Ref. [36] . 




Fig. 1. Density dependence of the strong part of the neutron- 
proton mass difference. The abscissa represents A = p/po, while 
the ordinate represents zlmnp*"^""^^ in units of MeV. 




Fig. 2. Density dependence of the electromagnetic part of the 
neutron-proton mass difference. The abscissa represents A = 
p/po, while the ordinate represents Aml'j^^'' in units of MeV. 



7 Results and discussion 

The medium-dependent effective masses and magnetic mo- 
ments of the nucleons are presented in Tab. [T] for a couple 
of values of the nuclear density. As expected, the effec- 
tive masses decrease with increasing density of the nu- 
clear medium. The most dominant contribution to these 
changes comes from the explicit medium (1 — Xp) factor 
in the expression of the mass functional M^p of the non- 
perturbed system 

Of course, it is well established that the effective mass 
of the nucleon decreases in the medium |571l38j . However, 
the in- medium behavior of the neutron-proton mass differ- 
ence is still under debate. Our results for the mass differ- 
ence due to the strong interaction, Z\mnp''^°"^\ and elec- 
tromagnetic interaction, Z\mn[f^'' , are presented in Fig. [1] 
and Fig. 2, respectively. 

^ In Ref. [21], where the Skyrme term is absent, this value is 
even equal to ~ 0.01 MeV. 



It can be seen that both parts, the strong and electro- 
magnetic one, have an almost linear dependence on the 
medium density and decrease with increasing the density. 
Note, however, that the absolute values of the changes are 
very small: i.e., the strong part of Amlp is almost flat 
(see Fig.[T]), while the change in the electromagnetic part 
is slightly more pronounced (see Fig.[2|). The difference in 

their behavior follows from the fact that zirnnp''""^'' has 
no explicit dependence on the medium functionals (see 

Eqs. ([IZl) and ((271)), whereas Am*Jp'^^ §B explicitly de- 
pends on the P-wave medium functional via the isovec- 
tor magnetic form factor ([M)) . However, there is a chance 
that the corresponding shifts may be more pronounced 
in isospin-asymmetric matter, due to the pertinent in- 
medium functionals in the isospin breaking term of the 
Lagrangian AC* (see Eq. dTT]) ) and the additional depen- 
dence of the medium functionals on the difference of the 
neutron-proton (distribution) densities of the surrounding 
environment, Sp — pn — Pp- 
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the decrease is about 0.85 MeV. One can conclude that in 
the Skyrme model, contrary to the result of the work [7j , 
Z\m*p remains positive - even at high densities - if the 
surrounding environment has an isospin-symmetric struc- 
ture. Note that in isospin-asymmetric matter this picture 
may change. 

The density dependence of the neutron-proton mass 
difference is approximately linear and can be parameter- 
ized as 

Z\m*p ~ 1.3 MeV - Cp/po- 

Our calculation shows that C ~ 0.05 MeV, which is almost 
negligible in comparison with the calculation \^ in the 
framework of QCD smu rules, C ~ 1.1 1.7 MeV. 



Fig. 3. Density dependence of the total neutron-proton mass 
difference. The abscissa represents A = p/po, while the ordinate 
represents Am^p in units of MeV. 
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Fig. 4. Density dependence of the following neutron- 
proton mass differences, normalized to their free space val- 
ues: Z\m*[,^^'/^"^ir\ Amli'''°''^^ / Amif°''^\ AmlJ Am^^ 
plotted by the dash-dotted, dashed and solid curves, respec- 
tively. The abscissa represents A = pI Pa, while the ordinate 
represents the three dimensionless ratios. 

In Fig.[3]the total neutron-proton mass difference Am*^^ 
in the nuclear medium is shown as function of the nu- 
clear density. For convenience, we present also, normalized 
to their free space values, the neutron-proton mass dif- 
ferences Z\mnif ^V^™n™\ Z\TOnp /Z\mnp™"^\ and 
Z\m*p/Z\mnp, in Fig. H] by the dash-dotted, dashed and 
solid curves, respectively. 

From Fig.[4]one can see that the change in the neutron- 
proton mass difference is only ~ 4% at normal nuclear 
matter density and that the total value Am\^ decreases in 
the nuclear medium as function of the density. This result 
is in qualitative agreement with the result of the quark- 
meson coupling model [7]. At the quantitative level, how- 
ever, they differ: the absolute value of the quantity Z\m*p, 
within the present approach, decreases about 0.05 MeV at 
normal nuclear matter density (see Fig.[3l), while in Ref. [7] 



8 Summary and outlook 

We have investigated the isospin-breaking effects for nu- 
cleons embedded into an isospin-symmetric nuclear envi- 
ronment. In order to calculate these effects, a medium- 
modified version of the Skyrme model, which also takes 
into account explicit isospin breaking, has been proposed. 
Our calculations within this framework show that, with 
increasing density, the total neutron-proton mass differ- 
ence as well as its strong and its electromagnetic part, 
separately, decrease only by a very small amount. 

A generalization of the present approach to finite nu- 
clei, i. e. to calculations of the type presented in Ref. [2U] , 
would allow to evaluate the isospin-breaking effects in mir- 
ror nuclei with the following qualifications: 

- the strong part of the neutron-proton mass difference 
would separately depend on the S- and P-wave medium 
functionals x^p, Xs,p s^-nd, consequently, would strongly 
depend on the parameterization of the latter ones [22] ; 

- the medium functionals would explicitly depend on the 
neutron- proton density distributions X^p" = X^p* {Pn + 
Pp; Pn — Pp) t22 , such that the main property of mirror 
nuclei 

could be taken into account in a natural way; 

- additional effects could arise due to an w-dependence 
of the polarization operators 77^'° = n^'°{u!,k) in 
Eq. ([TlPl . e.g. from the Weinberg- Tomozawa [391140] 
and the so-called range terms [U], see also Refs. [HJ 
S3]; 

- the evaluation of local isospin-breaking effects, accord- 
ing to the nucleon-position R as measured from the 
center of the nucleus [20] , would lead to additional pos- 
sibilities. 

The consequences of these additional modifications are 
the subject of future studies. 

Note that in the present work, the polarization operator, 
because of the form of Eq. ([6]) and the isospin symmetry of the 
surrounding environment, is simply given by the static formula 
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